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ABSTRACT 

We present a large set of merger simulations of early-type disc galaxies with mass 
ratios of 1:1 and 3:1 and 10% of the total disc mass in gas. The internal orbital 
structure and the kinematic and photometric properties of the remnants are analysed 
in detail and compared to pure stellar mergers. In contrast to the collisionless case 
equal-mass mergers with gas do not result in very boxy remnants which is caused by 
the suppression of box orbits and the change of the projected shape of minor-axis 
tube orbits in the more axisymmetric remnants. The isophotal shape of 3:1 remnants 
and the global kinematic properties of 1:1 and 3:1 remnants are only weakly affected 
by the presence of gas. 1:1 remnants are slowly rotating whereas 3:1 remnants are 
fast rotating and discy. The shape of the stellar LOSVD is strongly influenced by gas. 
Within the effective radius the LOSVDs of collisionless remnants have broad leading 
wings while their gaseous counterparts show steep leading wings, more consistent 
with observations of elliptical galaxies. We show that this change is also caused by the 
suppressed populating of box orbits and it is amplified by the formation of extended 
gas discs in the merger remnants which might eventually turn into stars. If elliptical 
galaxies have formed from mergers our results indicate that massive, slowly rotating 
boxy elliptical galaxies can not have formed from dissipative mergers of discs. Pure 
stellar (dry) mergers are the more likely candidates. On the other hand lower mass, 
fast rotating and discy ellipticals can have formed from dissipative (wet) mergers of 
early-type discs. So far, only unequal-mass disc mergers with gas can successfully 
explain their observed substructure. This is consistent with the revised morphological 
classification scheme of increasing importance of gas dissipation when moving from 
boxy ellipticals to discy ellipticals and then to spiral galaxies, proposed by Kormendy 
& Bender. 

Key words: methods: analytical - methods: N-body simulations - galaxies: elliptical 
and lenticular, cD - galaxies: formation - galaxies: evolution - galaxies: fundamental 
parameters 



1 INTRODUCTION 

Inspired by the work of iToomre fc Toomrel il972f) a large 
number of simulations have been performed to investigate 
whether giant elliptical galaxies can be formed by binary 
mergers of disc galaxies. If this formation is correct the 
properties of the merger remnants would have to be in 
agreement with all observed properties of elliptical galax- 
ies: their global kinematic and photometric properties as 
well as their fine-structure. Vi olent relaxation (iLvnden-Belll 
ll967l:lArad fc Johanssor]|2005D of the stellar and dark mat- 
ter components is the dominant dynamical process in the 
rapidly varying potential of an ongoing merger. However, it 
does not completely erase the information about the pro- 
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genitor galaxies. Therefore we should observe signposts of 
the formation process in the remnants, e.g. the morphology 
of the progenitor galaxies or the encounter geometry, both 
in the dynamics and the photometric properties. Elliptical 
galaxies, if they have formed this way, should show simi- 
lar properties, providing information about the properties 
of their progenitors. 

Collisionless simulations of equal mass merg- 



ers of disc galaxies have been studied in great de- 


tail by several authors (e.g. 


iGerhardl Il98ll: iBarnes! 


Il992t iHernauisd Il992l: iNaab fc 


Burkerd I2003L hereafter 


NB03 Gonzalez-Garci'a & Balcells 


120051 Ijesseit et alJ 


120051 hereafter JNB05. 


Naab & Truiillol 120061. 



iGonzalez-Garcia fc Balcells] boOa) ! The global proper- 
ties of the remnants are in agreement with observations 
of giant elliptical galaxies in the intermediate mass range, 
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e.g. equal mass remnants are slowly rotating, anisotropic, 
have boxy or discy isophotes. Unequal mass mergers are 
more isotropic and have discy isophotes (NB03). Merger 
remnants in general have phase space densities and surface 
density profiles that resemble obser ved ellipticals if bulges 
are added to the pro genitor discs ijHernauist et alJ 1199.4 
iNaab fc TruMc1l2006h . 

As soon as the kinematics of simulated merger remnants 
was investigated in more detail, an interesting disagreement 
with observed elliptical galaxies was revealed. The line-of- 
sight velocity distributions (LOSVD) within the effective 
radius of merger remnants in general show small asymmet- 
ric deviations from Gaussian shape. They tend to have a 
steep t railing wing (iBendo fc B arnes 200ft iNaab fc Burkertl 
l2001bT) . In contrast, most observed rotating elli pticals clearly 
show a steep leading wing in their LOSVDs teender et alJ 
11994) . Theoretically, axisymmetric, rota ting one-component 
syste ms in fact show such a behaviour iDehnen fc Gerhard! 
1994). This would indicate that ellipticals are very simple 
one component systems that did not form by mergers, It 
is, however, unlikely that elliptical galaxi es are such sim- 
ple systems (see e.g. lEmsellem et al.ll2004T) . An alternative 
explanation, based on photometric and kinematical obser- 
vations, is that rotat ing ellipticals contain e mbedded large 
scale stellar disc s fe.g. lRix fc Wmtell99CUScorza et alll998t 
iRix et al"l Il999l) . A superposition of two distinct compo- 
nents, e.g. a hot spheroidal bulge and a rotationally sup- 
ported cold disc, can also result i n a steep leading wing o f 
the LOSVD teender et alJll994T l. IrTaab fc Burkertl feOOlbD 
showed that embedding an exponential stellar disc artifi- 
cially in collisionless merger remnants would change the 
asymmetries in the LOSVDs, leading to a good agreement 
with observations. Those discs need a scale length similar to 
the effective radius of the bulge and 10% -20% of its mass. 

Such an extended disc naturally forms in merger rem- 
nants from high angular momentum gas located in the 
outer regions of the progenitors after it has been expelled 
along the tidal arms jNaab fc BurkerlJl2001al: lBarnesll2002t 
ISpringel fc Hernauisdl2005l: IWetzstein et al.ll2005ir Gas ini- 
tially located at smaller radii shocks and falls to the cen- 
ter of the remnant probably causing a starburst and/or 
feeding a super-massive black hole jMihos fc Hernauistl 
" 1996| : ISpringel et al-lliool . This gas does not contribute 



to the extended disc. It has been shown by iBekkil jl99i 
using simulations including gas and star formation that 
rotating lenticular galaxies with a disc like component 
can form from unequal-ma ss disc mergers. Furthermore, 
ISpringel fc Hernauistl (l2005t) has demonstrated with a simu- 
lation including star formation and feedback that a remnant 
with a dominant stellar discs can form self-consistently in an 
equal-mass merger of two gaseous discs. However, those ini- 
tial conditions might be more typical for high redshift discs 
and for such an extreme case starformation can not be ne- 
glected. 

In this paper we present simulations of disc galaxy 
mergers with small gas fractions of 10% which are more 
similar to low redshift discs. Our aim is to understand the 
influence of a small dissipative component (which will al- 
way exist, even if star formation is considered) on the global 
structure of mergers remnants. This will help to constrain 
the importance of additional physical processes li ke star for- 
mation, black-hole accretion and feedback (e.g. ICox et all 



l2005t ISpringel fc Hernauis^l2005l: iRobertson etai]|2006l) . It 
has been shown bv TBarnes^^fern^urstT ~ i^9rl) that the 
presence of gas influences the stellar component of the 
merging galaxies. Gas that accumulates at the center of 
merger remnants changes the central pote ntial and can in- 
crease the velocity dispersion of the stars. iRobertson et alJ 
(2006) have demonstrated this effect in merger simula- 
tions including star formation. In addition, a steeper cen- 
tral potential well results in a more axi symmetric cen- 
tral shape of the remnants teamed Il99ll . At the same 
time the fraction of stars on box orbits is significantly 
reduced and tubes becom e the dominant orbit family 
tearnes fc Hernauisdll996h . The most reasonable explana- 
tion for this behaviour is that systems with steep cusps 
in their po tential can not sustain a large population of 
box orbits llGerhard fc Binnevl Il985l ISchwarzschildl Il99: " 



[ Merritt fc Fridmanl 119961: IValluri fc Merritd Il998t iBarnei 

In this paper we present a new and more de- 
tailed analysis. With respect to resolution, statistical 
completeness and comparison to observations our study 
goes beyond previous investigations of similar simu- 
lations (e.g. iNegroponte fc White! Il983l iBamesI Il992l : 
iBarnes fc Hernauistlll996| : IBendo fc Barnesll2000fl . We show 
quantitatively how a dissipative component changes the 
LOSVDs and the photometric properties of the merger rem- 
nants and explain in detail how these changes are connected 
to the changes of the orbital content of the remnants (see 
JNB05 for collisionless remnants). In particular, we demon- 
strate that a dissipative component can explain the ori- 
gin of the observed asymmetries in the LOSVDs of ellip- 
tical galaxies. This paper is accompanied by a detailed two- 
dime nsional kinematic a nalysis of the remnants presented 
here J Jesseit et alJl200rj) . 

The paper is organised as follows. Section [5] presents 
the details of the simulations. A short overview of observ- 
able properties of ongoing mergers is given in Section|3] The 
intrinsic and projected shapes of the remnants and their or- 
bital content are discussed in Section 0] In Section |S] we 
show the projected photometric and kinematic properties 
of the remnants followed by a detailed investigation of the 
LOSVDs in Section A summary and conclusions follow in 
Section □ 



2 THE MERGER MODELS 

Disc galaxies were constructe d in dynamical eq uilibrium us- 
ing the method described bv lHernauistl J1993T) . The system 
of units was: gravitational constant G=l, exponential scale 
length of the larger progenitor disc = 1 (the scale height 
was h z = 0.2) and mass of the larger disc Md = 1. The 
discs were exponential with an additional spherical, non- 
rot at ingbjikje_mHT_jnass Mi, — 1/3, a Hernquist density 
profile llHernauisdll990D and a scale length rt = 0.2hd, and 
a pseudo-isothermal halo with a mass Md = 5.8, cut-off ra- 
dius r c = 10hd and core radius 7 = lh. The parameters 
for the individual components were the same as for the col- 
lisionless mergers presented in NB03. Scaled to the Milky 
Way, one length unit is S.bkpc, one velocity unit is 262km/ s 
and the mass unit is 5.6 x 10 10 Mq. For this study we re- 
simulated the full set of 1:1 and 3:1 mergers with an addi- 
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Figure 1. Time evolution of the gas mass inside a radius of 
r = 0.1 of the galaxies for 1:1 mergers (top panel) and the large 
(middle panel) and small (bottom panel) progenitor galaxies of 
3:1 mergers. The central region of each galaxy has been followed 
separately. As the galaxies merge the enclosed masses become 
identical. The total initial gas mass of the more massive discs is 



(see e.g. lNaab fc Burkerdl2001al : iBarnedbool iLi et alll2005l 

and references therein). 

We followed mergers of discs with mass ratios of r\ — 1 
and rj = 3 where r\ is the mass of the more massive galaxy 
divided by the mass of the merger partner. The equal-mass 
mergers were calculated adopting in total 440000 particles 
with each galaxy consisting of 20000 bulge particles, 60000 
stellar disc particles, 20000 SPH particles representing the 
gas component in the disc, and 120000 halo particles. We 
decided to use twice as many halo particles than disc par- 
ticles to red uce heating and instability effects in the disc 
components jNaab et al.lll999T) by encounters between halo 
and disc particles. For the mergers with r\ = 3 the param- 
eters of the more massive galaxy were as described above. 
The low-mass companion contained a fraction of 1/rj the 
mass and the number of particles in each component, with a 
disc scale-length (stars and ga s) of h = \/l/r/ ] as ex pected 
from the Tully-Fisher relation JPierce fc Tullvlll992l) . 

The N-body/SPH simulations were performed using the 
hybrid N-body/SPH treecode VINE (Wetzstein et al., in 
preparation) with individual timesteps. The gravitational 
forces were softened with a Spline kernel of h siav = 0.05. 
The minimal size of the Spline kernel used for computing the 
SPH properties, /isph, was fixed to the same value. Implic- 
itly, this procedure suppressed gas collapse on scales smaller 
than the softening scal e and prevents numerical instabilities 
feate fc Burkerj|l997f) . All simulation have been run on a 
cluster of 64 1.5GHz Sun CPUs at the Institute of Astron- 
omy in Cambridge. 

The initial discs were run in isolation for two dynami- 
cal times to allow the systems to finally settle into an equi- 
librium state. In the merger the galaxies approached each 
other on nearly parabolic orbits with an initial separation of 
30 length units and a pericenter distance of 2 length units. 
A study of orbits of merging da rk matter halos in cosmolog- 
ical large scale simulations by iKhochfar fc Burkerd j2006l) 
has shown that a significant number of the merging halos 
are indeed on parabolic orbits with a broad distribution of 
pericenter distances. In this study we focus on the usually in 
simulations considered regime of small pericenter distances. 
In selecting unbiased initial parameters for t he disc i nclina - 
tions we followed the procedure described bv lBarnesl lll99Sl) . 
The initial orientations for the discs were the same as in 
NB03, Table 1. The merger remnants were allowed to settle 
into dynamical equilibrium for approximately 30 dynami- 
cal time-scales after the merger was complete. Then their 
equilibrium state was analysed. 



tional gas component in the disc. We replaced 10% of the 
stellar disc by gas with the same scale length and an ini- 
tial scale height of h Ztgas = 0.1h z . The gas was represented 
by SPH particles adopting an isothermal equation of state, 
P = cjp, with a fixed sound speed of c B = 0.039 in velocity 
units, corresponding to c s w 10km /s if scaled to a Milky 
Way type galaxy. Assuming an isothermal equation of state 
implies that additional heat created in shocks, by adiabatic 
compression and feedback processes is radiated away im- 
mediately. It also implies that substantial heating processes 
prevent the gas from cooling below its effective tempera- 
ture, e.g. sound speed. Recently, a number of simulations 
have shown that using an isothermal equation of state is a 
reasonably good approximation to the ISM in disc galaxies 



3 COMMENTS ON MERGER DYNAMICS 
AND GAS INFLOW 

The properties of the initial disc galaxies change dramat- 
ically during the interaction. Tidal forces lead to the for- 
mation of tidal arms and trigger gas inflow to the cen- 
ter. The stellar systems are dynamically heated and finally 
merge into a new type of galaxy. Most investigations have 
focused on the properties of the final merger remnants. How- 
ever, recent high resolution observations of nearby inter- 
acting galaxies have made it possible to compare ongoing 
mergers directly to d ifferent phases of simulated interactions 
JPasvra et alJl2006|) . 
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It has not been investigated in detail during which 
phases of the mergers e.g. the central velocity dispersion 
or the effective radius adjusts to its final value. A detailed 
knowledge of how important observables evolve helps to con- 
strain m easurements of e.g. the mass-ratios of observed disc 
mergers dPasvra et ablEoOr?) . their gas accretion history, or 
their starformation rates. 

We measured the observables, in particular the effec- 
tive velocity dispersion of the progenitor discs, during the 
interaction by following every individual galaxy analysing 
snapshots of the merger in the orbital plane every 10 unit 
times (approximately every 10 half-mass rotation periods of 
the more massive disc). To perform th e analysis as consis- 
tent as possible with observations fsee lDasvra et alll2006l) 
we computed the effective radius r e s of every galaxy as 
the projected spherical half-mass radius of the stellar par- 
ticles within 5 length scales only taking particles of the 
galaxy itself into account. Thereby we avoided unrealisti- 
cally large values for r e s when the galaxies overlap. The ef- 
fective central stellar velocity dispersion <r e ff for each galaxy 
was then computed within 0.5r e ff taking all stellar particles 
into account. To follow the gas accretion onto the center we 
computed the total gas mass within a radius of r — 0.1, 
M sas (r < 0.1), which is two times the gravitational Spline 
softening length of the simulations, ft e = 0.05. 

In Fig. we show the amount of gas accreted onto the 
center of each galaxy. Every line corresponds to one disc 
galaxy. As soon as the galaxies merge (at t ps 70) the two 
lines of the progenitor galaxies join. In general equal mass 
galaxies merge faster than unequal mass galaxies. For 1:1 
mergers every galaxy contains M gas = 0.1. After the first 
encounter between 10% and 60% of the available gas in the 
progenitor discs is funneled to the center. In some cases ad- 
ditional gas, expelled form the partner galaxy, is captured 
already at early phases of the merger. The exact numbers 
vary with the encounter geometry. In general, gas transport 
is more effective if the spi n and the orbi tal angular momen- 
tum of the disc is aligned <BarneJl20o3) . After the galaxies 
have merged between 50% and 85% for 1:1 mergers (between 
40% and 75% for 3:1 mergers) of the gas has accumulated 
at the center. 

To test the effect of halo resolution on the gas inflow 
we have re-run two 1:1 and 3:1 mergers with very different 
geometries, respectively, changing the halo properties in two 
ways: First we have changed the softening of the halo to 
thaio = 0.09 and the gas softening to e gas = 0.03 to scale as 
the square root of the particle mass. This guarantees that 
the maximum gravitation al force exerte d from a particle is 
independent of its mass dDehnenl 120011) . The effect of this 
change on the gas properties was very small and we could 
not find a general trend. 

Secondly we used the original halo softening (ehaio = 
e — 0.05) and increased the mass resolution of the more 
massive halo to 348000 particles to get the same individual 
particle masses in all collisionless components. For this set of 
simulations we found that dissipational effects are slightly 
enhanced. The amount of gas driven to the center is in- 
creased by 5% to 15%. In Fig. we show the cumulative gas 
mass distribution of one 1:1 and on 3:1 merger in compar- 
ison to the original simulation. Most likely the higher halo 
resolution reduces the effect of artificial two-body heating 
by encounters of halo and gas particles. The global effect 
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Figure 3. Gas mass inside a radius of r = 0.1 of the galaxies for 
1:1 mergers (top panel) and the large (middle panel) and small 
(bottom panel) progenitor galaxies of 3:1 mergers as a function 
of nuclear separation. 



on the gas distribution is, however, small and the stellar 
remnants are only weakly affected (see Section 0J. Based on 
these results and to be fully consistent with the already per- 
formed collisionless simulations we have used a gravitational 
softening of e = 0.05 for all components. 

In Fig. |2]we show the final gas distribution for 1:1 and 
3:1 mergers with three different encounter geometries. For 
1:1 remnants the gas can be distributed in warped discs 
or more irregular structures like polar rings. 3:1 remnants 
show a more regular disc structure. The formation of ex- 
te nded gas discs in 3:1 merg er remnants has been pr e dicted 
bv lNaab fc Burkertl feOOlal) and we refer to iBarneJ i2002l) 
who has studied the dynamics of disc formation extensively. 
In Fig. [3] we show the central gas content as a function of 
nuclear separation of the merging galaxies. The main ac- 
cretion phases are after the first encounter and at the fina l 
merger of the galaxies (see e.g. iMihos fc Hernauistlll996T) . 



The influence of gas on the structure of merger remnants 5 





Figure 2. Final gas distribution for 1:1 and 3:1 remnants with three different geometries A, B, and C seen along the principal axes of 
the main stellar bodies. The box length is 8 unit lengths in vertical direction and 18 unit lengths in horizontal direction. Most remnants 
have embedded gas discs with a more regular disc structure for 3:1 remnants. 



It is important to note for comparisons with observations 
that the central gas content differs significantly before and 
after the first encounter although the galaxies can have the 
same nuclear separation. This fact could be used as addi- 
tional diagnostics to predict in which state of the merger 
observed galaxies are. The effective velocity dispersion ver- 
sus nuclear separation is plotted in Fig. 2] Similar to the gas 
content the dispersion rises in two phases: after the first en- 



counter and, to its new equilibrium value, in the final merger 
phase. In combination with the evolution of the effective 
radii these results have been used to test dynamical mass 
estimates of observ ed nearby ULIRGs in interacting pairs 
JPasvra et alj f2006) which all appear to have mass ratios 
between 1:1 and 3:1. Our results also indicate that merging 
disc galaxies might already fall on the observed black-hole 
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Figure 6. Axis ratios of the stellar remnants for the simulations with and without gas for mass ratios 1:1 (black dots) and 3:1 (red 
diamonds) Left panel: ratio between the short and the long axis (c/a). Middle panel: ratio between the intermediate and the long axis 
(b/a). Right panel: triaxiality of the stellar remnants. Stellar remnants with a gas component are less triaxial with a similar maximum 
flattening (c/a). 



mass - sigma relation jTremaine et alJl2002l) soon after their 
first encounter (Dasyra et al., submitted). 



4 INTRINSIC SHAPES AND STELLAR 

ORBITS OF THE MERGER REMNANTS 

The intrinsic shape of a mass distribution is defined by the 
ratio of its three principal axes. They were determined by 
diagonalising the moment of inertia tensor of each merger 
remnant. The particles were binned according to their bind- 
ing energy. That ensures that the subsets of particles follow 
the s tructure of the remnant naturally jWeil fc Hernauistl 
1996). The triaxiality parameter T is defined as 



where a, b, and c are the long, intermediate and minor-axis, 
respectively. The shape of merger remnants is closely re- 
lated to their intrinsic orbital structure as shown by JNB05. 
In general minor-axis tubes are more heavily populated in 
oblate remnants and box orbits are most abundant in rem- 
nants triaxial remnants with T = 0.5 (see Fig. 5 and Fig. 6 
of JNB05). In FigEJwe show how the presence of gas influ- 
ences the intrinsic shapes of the stellar components of the 
merger remnants. The triaxiality is lowered for almost ev- 
ery remnant due to the influence of gas. 3:1 remnants with 
gas only reach a maximum triaxiality of T = 0.2, while in 
the collisionless case they can have triaxialities as large as 
T — 0.5. Although the triaxiality of the 1:1 remnants is also 
lowered substantially the effect is smaller than for 3:1 rem- 
nants. The more violent merging process can to some extent 
counter the dissipational influence of the gas. Gas rather 
makes the remnant more axisymmetric than more flattened 
(Fig|HJ|. The ratio between the short and the long axis of 
the triaxial body (c/a) is only increased by ~ 10%, making 
the remnant slightly more spherical. In contrast, for almost 
all remnants b/a is now close to unity, i.e. the systems are 
axisymmetric. 

The classification of orbits follows the procedure pre- 
sented in JNB05 (and re ferences therein) and i s base d on 
the method developed bv ICarpintero fc Aguilarl jl998l) . We 
will repeat only the most important steps. After the merger 



remnant had settled into equilibrium we froze the particle 
distribution and computed the potential using the Self C on- 
sistent Field (SCF) method Jrlernciuist fc Ostrikerll992l) . As 
initial condition for the integration of the orbits we took the 
positions and velocities of the particles at the final snapshot 
of the simulation. With this method all orbit classes which 
exist in general triaxial potentials can be identified. The 
orbits were classified as minor-axis tubes, outer major-axis 
tubes, inner major-axis tubes, boxes and boxlets. About 5% 
-10% of all orbits in every remnant could not be classified 
with this procedure. 

Every orbit class has distinct kinematical properties and 
shapes. Minor-axis tubes are the backbone of oblate and 
disc-like systems and have a non-vanishing angular momen- 
tum around the short axis of the potential. Major-axis tubes 
are found in prolate or nearly spherical systems with a non- 
vanishing angular momentum around the long axis. Boxes 
and boxlets have no mean angular momentum and can be 
found close to the center of the potential. 

We compared the population of each orbit class aver- 
aged over all 1:1 and 3:1 mergers, respectively, for all merger 
remnants with and without gas. The result is shown in Fig. 
[7| In both, 1:1 and 3:1 mergers, gas drastically reduces the 
fraction of box orbits while the fraction of minor-axis tubes 
increases by a factor of 2-3 (the slightly larger central gas 
fractions of the simulations with higher halo resolution dis- 
cussed in the previous section results in an increase of the 
tube orbit fraction of 3-5%). This confirms the observed 
trends towards an axisymmetric shape. Major-axis tubes 
are also de-populated which is consistent with the reduc- 
tion of the prolateness. Box orbits are still found at the 
smallest radii, but minor-axis tubes in general start to be 
the dominant orbit class at smaller radii than in the dissi- 
pationless mergers ( FigEl . This finding i s in ag reement with 
iBarnes fc Hernauisd <ll996l) and lBarnesI <ll998t) . 

The superposition of all stars, which are on different 
orbits, determine the projected photometric and kinematic 
properties of the remnants which are investigated in the fol- 
lowing sections. 
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Figure 9. Typical projected surface density of stellar particles on minor-axis tube orbits for a triaxial collisionless 1:1 remnant (upper 
row) and its more oblate counterpart with gas (bottom row). From left to right we show projections along the short, intermediate, and 
long axis. The box length is PS 2r e ff . Tube orbits, which dominate around the effective radius, support a boxy/peanut isophotal shape in 
the triaxial potential of the collisionless remnant. In the more axisymmetric potential of the remnant with gas the shape is less boxy or 
even discy. 



5 GLOBAL PHOTOMETRIC AND 
KINEMATIC PROPERTIES 

We have performed the isophotal and kinematic analysis 
of 500 random projections of every remnant following the 
procedure presented in NB03. An artificial image of each 
projected remnant was created by binning the stars within 
the central 10 length units into 128 x 128 pixels, smoothed 
with a Gaussian filter of standard deviation 1.5 pixels. The 
isophotes and their deviations from perfect ellipses were then 
determined using a data reduction package kindly provided 
by Ralf Bender. In this section the analysis was only per- 
formed for the stellar particles. In the following we refer 
to the stellar particle distribution as collisionless/without 
gas or dissipative/with gas depending on whether gas was 
present or not. 

The characteristic ellipticity e e ff for each projection is 
denned as the isophotal ellipticity at 1.5r e . In Fig. 1101 wc 
show the normalized histograms of the ellipticity distribu- 
tion for the stellar 1:1 and 3:1 remnants with and without 
gas. The distributions look similar for 3:1 remnants, whereas 
1:1 remnants with gas have more projections with small el- 
lipticities as expected if the main stellar body is more ax- 
isymmetric. Still, the distributions peak at similar elliptici- 
ties as mainly the ratio of the intermediate to the long axis, 
b/a, is affected by the presence of gas (see Section 0. 

The effective a4-coefncient, a4 e g, was computed as the 
mean value of between 0.25r e and 1.0r e , with r e being the 
projected spherical half-mass radius. In contrast to NB03 we 
do not use the maximum value of 04 in case of a peaked dis- 
tribution which did, however, not change the results. The 
normalized histograms of the a4 c ff distribution for all 1:1 
and 3:1 remnants with and without gas is shown in Fig. 
1111 In contrast to the collisionless case the 1:1 remnants 



with gas do not shown a significant number of boxy projec- 
tions any more. Their deviations from perfect ellipses has 
decreased. For 3:1 remnants the effect of the gas is weaker. 
There are now more projections which are significantly discy, 
with a4 c ff > 2. 

The most important physical reason for the lack of boxy 
projections for mergers with gas is the different behaviour of 
minor-axis tube orbits in axisymmetric and triaxial poten- 
tials. As we have shown in Section 3] collisionless 1:1 rem- 
nants are more triaxial whereas stellar 1:1 remnants with 
gas are more axisymmetric. Minor-axis tube orbits in triax- 
ial potentials can support a boxy/peanut isophotal shape in 
the projection along the minor axis and, more importantly, 
along the long axis (see Fig. 11 in JNB05). Minor-axis tubes 
are the dominant orbit family around the effective radius 
of the collisionless remnants (see Fig. |HJ and therefore are 
largely responsible for the over all isophotal shape (see NB03 
and JNB05). In the more axisymmetric potentials of the 
remnants with gas the minor-axis tubes look more ellipti- 
cal/less boxy or even discy in all projections. This behaviour 
is qualitatively demonstrated in Fig.|§]for a typical 1:1 rem- 
nant. In addition, the fraction of box orbits and boxlets, 
which can support a boxy shape at the center of collision- 
less remnants (see Fig. 11 in JNB05) is significancy reduced 

(Fig. El. 

The central velocity dispersion cro of every remnant was 
determined as the average projected velocity dispersion of 
the luminous particles inside a projected galactocentric dis- 
tance of 0.2r e and we defined the characteristic rotational 
velocity along the major axis, « ma j, and the minor-axis, « m in, 
as the rotational velocity at 1.5r e and 0.5r e , respectively. 
The amount of minor-axis rotation was parameterized as 



« 1 / v^> 



un feinnevlll985T) . Minor-axis rotation in 
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Figure 12. Global photometric and kinematic properties of the stellar components of 1:1 and 3:1 merger remnants with and without 
gas. The contours indicate areas of 50% (dotted line), 70% (thin line) and 90% (thick line) probability to detect a merger remnant with 
gas with the given properties. For comparison, the shaded areas indicate the 50% (dark grey) and 90% (bright grey) probabilities for 
collisionless remnants. Data for observed boxy (filled squares) and discy (open diamond s) ellipticals hav e been kindly provided by Ralf 
Bender. Data for local merger remnants (open green squares) are taken from lRothberg fc Joseph! a2006f) . From top to bottom we show: 
cllipticity, e c ff, versus i> ma j/ou (the theoretical value for a spheroid flattened by rotation is shown by the dashed line); effective isophotal 
shape, a4 c £f, versus e cf f , versus the amount of minor-axis rotation, v m in/ v^niaj v 'min > an< ^ versus the anisotropy parameter,(i> ma j/(7o)*. 
The stellar components of remnants with gas (in particular 1:1) do not show boxy isophotes any more and the 3:1 remnants with gas 
are more discy and show less minor-axis rotation. 
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Figure 4. Effective stellar velocity dispersion, <r e ff , of the galaxies 
for 1:1 mergers (top panel) and the large (middle panel) and small 
(bottom panel) progenitor galaxies of 3:1 mergers as a function 
of nuclear separation. 



elliptical galaxies, in addition to isophotal twist, has been 
suggested as a sig n for a triaxial shape of the main body of 
elliptical galaxies JWagner et alJll98St iFranx et alJll99ll) . 

In Fig. 1121 we summarize the properties of the stel- 
lar components of the 1:1 and 3:1 merger remnants with 
gas in comparison to the remnants without gas. The re- 
maining gas is not included in this analysis. In the up- 
per row we show the location of the remnants in the 
iWj/ffo — £eff plane. There are no fundamental differences 
between remnants with and without gas. The dichotomy 
found bv lBarnesl lll99Sl) . iNaab et alJ i 19991) and NB03 that 
1:1 mergers produce slowly rotating remnants and 3:1 rem- 
nants produce fast rotating remnants is not changed by the 
presence of gas. In general the properties of 1:1 and 3:1 
remnants are in very good agreement with recent obser- 
vations of (predominantly discy) nearby merger remnants 
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Figure 5. Comparison of the cumulative gas mass versus radius 
for a 1:1 (solid line) and a 3:1 merger (dashed line) and their 
counterparts with higher resolution in the halo component (dot- 
ted and dash-dotted line). The shaded area indicates the central 
part of the remnants. A better resolved halo results in a more con- 
centrated gas distribution. The central gas content can increase 
by 15%. 
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Figure 7. Mean orbit fractions for the 40% most tightly bound 
particles for all 1:1 (circles) and 3:1 (diamonds) remnants without 
(open symbols) and with gas (filled symbols). 



bv lRothberg fc Joseph! i2006T) and discy and boxy elliptical 
galaxies. 

The distribution of the remnants in the e e g — &a plane 
is shown in the second row. As shown above the stellar 
components of 1:1 remnants with gas do not appear boxy 
any more and can not explain the formation of boxy ellipti- 
cals. For individual simulated m ergers this trend was already 
fo und bvlBekki fc Shioval l)l997t) and thereafter rediscovered 
bv lSpringell (120001) . Without assuming additional physics for 
star formation all formerly boxy remnants are now close to 
elliptical (see Fig. I HP - The agreement of 3:1 remnants with 
nearby merger remnants and discy elliptical galaxies is sig- 
nificantly better. 

In the third row of Fig. 1121 we compare the amount of 
minor-axis rotation as a function of the isophotal shape of 
the collisionless remnants and the remnants with gas. There 
is a weak trend for remnants with gas to show less minor- 
axis rotation than the collisionless counterparts. This finding 
can be understood as the fraction of major-axis tubes, which 
are mainly responsible for minor-axis rotation, is on average 
reduced due to the influence of gas (see Fig. . 

The last row of Fig. 1121 shows the location of the rem- 



10 Thorsten Naab, Roland Jesseit & Andreas Burkert 



0.8 



fx 

g 0.4 



0.2 



Boxlet 
C>-f> Minor Axis Tubes 
00 Outer Major Axis Tubes 
<r<] Inner Major Axis Tubes 
& Boxes 




1:1 gas 




3:1 no gas 




3:1 gas 




0.12 



0.10 - 



1 


1 


no gas 


1 


1 


gas 


3 


1 


no gas 


3 


1 


gas 



c 
2 

Oh 

c 

o 
o 

(0 

Sh 



0.08 - 



0.06 - 



0.04 



0.02 



0.00 




0.0 0.2 0.4 0.6 



0J 



Figure 10. Distribution of the projected effective ellipticities, 
e e ff for the stellar body of 1:1 remnants with (solid) and with- 
out (dashed) gas and for 3:1 remnants with (dashed-dotted) and 
without (dotted) gas. 
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Figure 8. Fraction of particles on different orbits in radial bins 
in units of r e for a typical 1:1 remnant without and with gas (two 
upper panels) and a typical 3:1 remnant with and without gas 
(two lower panels). Remnants with gas are more dominated by 
tube orbits. 



Figure 11. Distribution of the projected effective isophotal 
shape, ci4 c g , of the stellar body for 1:1 remnants with (solid) and 
without (dashed) gas and for 3:1 remnants with (dashed-dotted) 
and without (dotted) gas. 
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Figure 13. Kincmatical properties of a typical 3:1 merger rem- 
nant (collisionless, stars, stars+gas) . The shaded area indicates 
the radial range where the local correlations between /13 and 
i>fl t /(jfi t have been measured (see Fig, 1141 , If gas is present during 
the merger, h% changes sign. 

nants in the (vmaj/oo)* — 04 plan e. Here (twi /oo) * is the 
tradit io nal anisotropy paramet er feinnevlll97Sl see iBimievI 
( 2005*): iBurkert fc N aab (2005) for a revised version of the 
anisotropy parameter and its application to N-body simula- 
tions). On average the 1:1 remnants with gas have slightly 
larger values of («maj/0"o)* • Boxy and anisotropic remnants 
which we found for collisionless remnants did not form. In 
addition, we find more projected gas-remnants with small 
(Vmaj/00)* but discy isophotes than for the pure collision- 
less remnants. In this regi me no observed elliptical galaxies 
jRothberg fc Josephl200rl can be found. H owever, there ar e 
some SO galaxies with similar properties iRix et al.lll992l) . 
A possible connection will be investigated in a subsequent 
paper. 3:1 remnants with gas have similar properties to their 
collisionless counterparts. Remnants of mergers with gas on 
average appear more discy and are in better agreement with 
observed discy ellipticals and nearby merger remnants. 
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local v/a 

Figure 14. Local correlation between I13 and v/a for the same 3:1 
remnants as in Fig. 1131 The collisionless remnant shows an anti- 
correlation between for low v/a and a correlation for the largest 
values. For the stellar remnant of the merger with gas I13 and 
v/a is anti-correlated. Including the gas in the analysis leads to 
a stronger anti-correlation in good agreement with observations. 
The range of the observational data is indicated by the shaded 
area. 



6 LOSVD ANALYSIS OF THE MERGER 
REMNANTS 

To measure the LOSVDs of a merger remnant we shifted 
the densest region of every two-dimensional projection to 
the origin and placed a slit with a width of 0.4 unit lengths 
along the apparent long axis of each projected remnant. The 
slit was then subdivided into grid cells of 0.15 unit lengths. 
Thereafter we binned all particles falling within each cell in 
velocity along the line-of-sight. The width of the velocity 
bins was set to a value of 0.2 for line-of-sight velocities v\ os 
in the range —4 ^ v\ os ^ 4. This resulted in 40 velocity 
bins over the whole velocity interval. Using the binned ve- 
locity data we constructed line-of-sight velocity profiles for 
each bin along the slit. Subsequently we parameterized devi- 
ations fro m the Gaussian shape using Gauss-Hermite basis 
functions (lGerhardll993l:lvan der Marel fc Franxll993r) . The 
kinematic parameters of each profile (erst, ffl t, ^3, hi) were 
then determined by least squares fitting (see ICretton et alJ 
l200ll) . The large number of simulated stellar particles (> 
100000) guaranteed that at least 1000 particles fall within 
each bin inside the effective radius. We have analysed three 
sets of particle distributions: all star particles of the colli- 
sionless simulation, the star particles of the simulation with 
gas and, assuming that all gas particles have transformed 
into stars, stars and gas particles of the simulation with gas. 
In the following we refer to these sets of particles as: colli- 
sionless, stars, and stars+gas. 

As a prototypical example Fig. ll3l shows the observables 
«flt, ""fit, ffit/ofit, an d hi together with the bootstrapping 
errors as a function of radius for a 3:1 merger remnant. We 
compare the kinematics of the collisionless remnant with 
the stars and stars+gas of the simulated remnant with gas. 
The line-of-sight velocity is only weakly affected by the pres- 
ence of gas. There is a slightly steeper gradient at the cen- 
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Figure 16. Upper row: Local correlation between /13 and v/o~ between 0.25 and lr e along the major axis for 1:1 and 3:1 merger remnant s 
with and without gas. The contours are as in Fig. 1121 The vertical lines indicate the location of observational data bv lBender et alJ 
The collisionless 1:1 remnants (left panel) rotate very slowly. The asymmetry of the LOSVD of the stars (/13) increases if gas is added 
to the simulations. Collisionless 3:1 remnants do not show the observed anti-correlation (right panel, shaded area). The corresponding 
stellar components of the mergers with gas (contours) have larger maximum values of v/o~ and show stronger asymmetries with a clear 
indication for a change in the tilt. Second row. Including the gas (assuming it has transformed into stars after the merger was complete) 
the distribution for equal mass remnants becomes weakly tilted (left panel) whereas the distribution for 3:1 remnants (right panel) shows 
a clear tilt and extends to a local v/a ss 1 in good agreement with observations. Third row: Same as above but for /14. Most collisionless 
remnants (shaded areas) have a positive local values of /14 whereas the remnants with gas show equal amount of positive and negative 
values (contours). For 3:1 remnants there is a trend for regions with higher v/cr to have negative /14. Bottom panel: Adding gas to the 
analysis results in an over-all shift to more negative values of /14. 
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Figure 15. Typical normalised velocity profiles at Kb 0.75r e ff of 
the dominant population of Z-tubes and box orbits for a colli- 
sionless 3:1 remnant (upper panel) and the stellar component of 
its counterpart with gas (lower panel). Z-tubes have a negative 
h%. In combination with the large fraction of box orbits of the 
collisionless remnant the overall LOSVD shows a broad leading 
wing. The remnant with gas is dominated by Z-tubes and has 
a steep leading wing, in agreement with observations of rotating 
elliptical galaxies. 



ter for the merger with gas. The velocity dispersion pro- 
file (stars+gas) shows a clear signature of the extended gas 
disc: it decreases faster with increasing radius. At the center 
where the gas has formed a cold dense blob the differences 
in the dispersion are most significant. 

The most interesting feature can be seen for /13 which 
characterises the asymmetry of the LOSVD. In the colli- 
sionless case /13 and vat have the same sign - and nearly 
the same slope - inside half an effective radius. Towards 
one effective radius the slope for /13 changes and at larger 
radii vat and /13 are anti-correlated (Fi g. I13| . This be- 
haviour is typical for all 3:1 remnants (iNaab fc Burkertl 
l2001bl : iBeiido fe Barnesll2000l) and it is in contradiction to 
observed elliptical galaxies feender et al.l ll994) where /13 
and ffit/ofit are anti-correlated inside the effective radius 
(indicated by the shaded area in Fig. [T]J. The slope of the 
correlation between /13 and vat I oat changes if gas is present 
during the merger. In the region of interest /13 and vat/oat 
are now weakly anti-correlated for the stars and strongly 
anti-correlated for stars+gas (see Fig. 1141 . In Fig. 1151 wc 
show the velocity distribution and the Gauss-Hermite fits 
for box orbits and Z-tubes, which dominate the remnants. 
The collisionless 3:1 remnant is compared to its counterpart 
with gas (same initial orientation). For both remnants the 
Z-tubes have a negative /13 which corresponds to a steep 
leading wing. In combination with the large fraction of box 



orbits the resulting L OSVD of the collis i onless remnant has 
a broad leading wing. iBendo fc Barnesl i200Cj) have shown 
a very similar distribution for their remnants. The stellar 
component of the remnant with gas has a much lower frac- 
tion of box orbits and the resulting LOSVD is dominated 
by the steep leading wing of the minor-axis tubes. 

The fourth order coefficient, hi, characterizes a LOSVD 
that is more peaked than a Gaussian if positive and less 
peaked if negative (Fig. 1131 . For the 3:1 remnant chosen 
here hi appears to be positive at the center and negative 
around the effective radius, a feature that is less prominent 
for simulations with gas. In this paper we show the results 
for hi for completeness, however we will not discuss it in 
greater detail. 

In Fig. 1161 we have summarised the local correlations 
between /13 and v/a in the range of 0.25 to 1 r e along the 
projected major axis for 500 random projections for all 1:1 
and 3:1 merger remnants with and without gas. The colli- 
sionless 1:1 remnants (left panel) rotate very slowly and are 
consistent with observations. The asymmetry of the LOSVD 
of the stars (/13) increases if gas is added to the simulations. 
Rotating collisionless 3:1 remnants do not show the observed 
anti-correlation (right panel, shaded area). The correspond- 
ing stellar components of the mergers with gas (contours) 
have larger maximum values ofv/o and show stronger asym- 
metries with a clear indication for a change of tilt. Includ- 
ing the gas (assuming it has transformed into stars after the 
merger was complete), the distribution for equal-mass rem- 
nants becomes weakly tilted whereas the distribution for 3:1 
remnants shows a clear tilt and extends to a local v/a ~ 1 
in good agreement with observations. 

The lower panels of Figure Hol show the correlation with 
hi. Most collisionless remnants have positive local values of 
hi whereas the remnants with gas show equal amount of 
positive and negative values. For 3:1 collisionless remnants 
there is a trend for regions with higher v/a to have negative 
hi. Adding gas to the analysis results in an over-all shift to 
more negative values of hi. 



7 SUMMARY AND DISCUSSION 

We have presented a statistical sample of simulations of disc 
galaxy mergers with a 10% fraction of gas in the progenitor 
discs. The properties of the merger remnants have been com- 
pared to the properties of a similar set of simulations without 
gas. The effect of star-formation was not included. The pres- 
ence of a dissipative component changes the shape and the 
orbital content of the stellar component of the merger rem- 
nants. The fraction of box and boxlet orbits which dominate 
the inner parts of collisionless remnants is significantly re- 
duced if gas is included. The fraction of outer major-axis 
tubes is reduced as well and inner major-axis tubes dis- 
appear completely. In remnants with gas most stars move 
on minor-axis tube orbits. The change of orbits is caused 
by gas that settles at the center of the remnants where it 
deepens the potential well at the same time making it more 
axisymmetric. In this environment box orbits cannot exist 
as they are only supported in triaxial potentials. We find 
that the intermediate axis of the stellar distribution is most 
strongly affected and the stellar remnants with gas are more 
axisymmetric. These results are in good qualitative agree- 
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ment with th e findings of iBarnes fc Hernauisd (119961) and 
iBarnej (llQQSTl . However, with respect to resolution, statisti- 
cal completeness and comparison to observations our study 
goes beyond previous investigations. 

The isophotal shape of equal-mass remnants is strongly 
affected by gas. Around the effective radius minor-axis tubes 
- as well as box-orbits and boxlets at smaller radii - in 
the triaxial potential of collisionless remnants support boxy 
isophotal shapes. In the remnants with gas the fraction of 
box orbits is reduced and the now dominant tube orbits ap- 
pear more elliptical or even discy. Statistically, equal-mass 
mergers with gas do not produce boxy remnants, but they 
still have a small anisotropy parameter which is in con- 
flict with observations of ellipticals. This is only valid in 
the limiting case of no starformation during the merger. 
Realistically, equal-mass disc mergers li ke nearby UL IRGs 
do ex perience bursts of star formation jGenzel et al.lll993 . 
l200lT) . It is, however, still unclear how much of the available 
gas is transformed into stars at which stage of the merger. 
For indivi dual equal-mass merge r si mulations i t has been 
shown by iBekki fc Shioval (119971) and ISpringell (|200CD that 
boxy remnants can be formed if the star-formation efficien- 
cies were chosen to be high lead ing to an evoluti on that 
is similar to the collisionless case llNaab et alJll999h . Disky 
remnants form for low efficiencies similar to the results with 
gas presented here. 3:1 rem nants with gas are slightly more 
discy (see e.g. [Bekki 1998 for an unequal-mass merger sim- 
ulation including star formation) but in general resemble 
their collisionless counterpart and are in go od agreement 
with observations of d iscy elliptical galaxies feender et alJ 
ll988|; |Hao_g_Lal] |2006j) as well as nearby merger remnants 
jRothberg fc Josephl2006l) . We expect that the effect of star 
formation on the isopho tal shape of unequal m ass remnants 
is much weaker (see e.g. iBournaud et al.ll20o3) . 

The shape of the LOSVD of the stars is significantly in- 
fluenced by the presence of gas. Collisionless remnants have 
LOSVDs which are close to Gaussian or have broad leading 
wings. In contrast the stellar component of mergers with 
gas shows are clear tendency for steep leading wings. For 
remnants with a significant amount of rotation (e.g. 3:1 rem- 
nants) this results in an anti-correlation between /13 and v/a, 
especially if the gas component that has formed extended 
discs is included in the analysis after the merger is com- 
plete. This is in good a greement with obse rvations of discy, 
fast rotating ellipticals feender et al]ll994r) and might indi- 
cate low star-formation efficiencies during the merger event 
itself. The LOSVDs of equal-mass merger remnants are also 
consistent with boxy, slowly rotating ellipticals. 

If mergers of discs galaxies have indeed formed ellip- 
tical galaxies in the not to recent past o ur re sults con- 
firm the scenario of iKormendv fc Bender! il996fl that gas 
dissipation becomes more important the more discy the 
isophotes of early-type galaxies are and the faster they 
rotate. Some equal-mass merger remnants with gas still 
have a small anisotropy parameter and discy isophotes. 
Either this conflict with observations can be solved by 
efficient star formation during the merger or it indi- 
cates that massive, boxy ellipticals formed in dissipation- 
less mergers from predominantly stellar progenitors which 
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cent high resolution direct numerical simulations of the for- 
mation of field elliptical galaxies from cosmological initial 
conditions indicate that early-type galaxies after they have 
formed at high redshift during a phase of intensive merging 
can thereafter gro w by accretion of mainly stellar satellites 
or minor mergers jNaab et alj|2005h . 

Combining the results from NB03 and the present study 
we can conclude that mergers of typical disc galaxies with 
bulges can have formed giant elliptical of low and interme- 
diate mass in the past. The agreement of kinematic and 
photometric properties of 1:1 and in particular 3:1 disc 
mergers with nearby merger remnants is striking. In com- 
bination with result s from the simulations presented here 
iDasvra et alJ (|2006) have been able to show that ongo- 
ing disc mergers with ULIRG activity have mass ratios 
between 1:1 and 3:1. Furthermore the properties of the 
simulated remnants are in good agreement with the kine- 
matics and isophotal sha pes of nearby merger remnants 
iRothber g fc Joseph) |20 06l) . A parti cularly interesting re- 
sult is that lRothbere&^osephl <l2006h find merger remnants 
which are discy and anisotropic, a regime that is populated 
by simulated remnants but not by virialised old elliptical 
galaxies (see their paper for a detailed discussion). The di- 
rect comparison of real nearby mergers with merger simu- 
lations are a powerful tool to constrain theories on ellipti- 
cal galaxy formation and the effect of st ar formation and 
black-hole formation jPasvra et alJ 120061) in the local uni- 
verse which then can be applied to simulations at all redshift 
ranges. 

In the present simulations the conversion of gas into 
stars has not been included and only progenitors with a 
small gas fraction have been considered as otherwise neglect- 
ing star formation would hardly be justified. We find impor- 
tant signatures of the small gas component in the merger 
remnants. The question of how stars form and how stellar 
energetic feedback or central black hole heating affect the 
multi- phase interstella r medium is still poorly understood 
(for a summary see e.g. lElmegreen fc Scalol20o1h . leading to 
large uncertainties in models of galaxy formation and evo- 
lution. Recent simulati on e.g. bv lSpringel et ail (l2005h and 
iDi Matteo et all <l2005l) however demonstrate their impor- 
tance in understanding the formation of elliptical galaxies 
in greater details. The orbital content of the remnants will 
be influenced by the timing of black-hole formation and the 
efficiency of feedback which determines the amount of gas 
that can make it to the center leading to the effects pre- 
sented here. It will be interesting in future work to explore 
the effect of these processes on the internal orbital struc- 
ture and projected properties of elliptical galaxies. From a 
kinematic point of view disc merger remnants appear very 
similar to observed ellipticals, although questions regarding 
the age and metalicity of the stellar populations have to 
be addre ssed in the future. Anal ytical models of disc forma- 
tion fe.g. lNaab fc OstrikeJ2006T) which contain all the infor- 
mation about stellar ages and metallic itites of the pro geni- 
tor d i scs and semi-analytical modeling iKhochfar fc Burkerd 
l2005t IKhochfar fc Silkll2006ft in comb ination with detailed 
merger simulations jNaab et alJl2006h can be used to place 
further constraints on the disc merger hypothesis. 
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